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We present results for the semi-leptonic and leptonic decays of B mesons. These non-perturbative matrix 
elements are important for constraining the CKM matrix. Results are presented for the pseudoscalar and vector 
decay constants, as well as flavour breaking ratios and heavy quark symmetry relations. We consider the chiral 
and momentum dependence of the semi-leptonic form factors of the decay B — > ttIu and the soft pion relation 
(SPR) on the lattice. These calculations were performed in the quenched approximation at two values of the 
coupling with non-perturbatively O(a) improved action and currents. 



1. INTRODUCTION 

The determination of decay constants and form 
factors arising from the decay of B and D mesons 
is of great importance to phenomenology, in par- 
ticular in constraining the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix. These quantities can 
be calculated on the lattice. 

The action used in this calculation is the 0{a) 
non-perturbatively improved action where 
csw and the coefficients needed to improve 
the currents, are determined non-perturbatively 
where possible. The calculation was done at 
two values of the coupling, j3 = 6.0 and j3 = 
6.2. The simulation parameters are shown in 
Table [|. We simulate with three light quark 
masses around strange and four propagating 
heavy quarks around charm. Heavy Quark Sym- 
metry (HQS) is used to motivate the form of the 
heavy extrapolation to the B scale. The light 
quark mass parameters (« C rit etc) have been de- 
termined by the UKQCD collaboration ||. The 

Table 1 

Simulation parameters. 





= 6.2 


13 = 6.0 


Volume 


24 a x 48 


16 3 x 48 


csw 


1.614 


1.769 


A^configs 


216 


305 


a-^ra) GeV 


2.913(10) 


2.119(7) 



vector and axial currents used in this calculation 



are improved and renormalised as follows: 
V* = Z v (l + b v am q )(v^+ac v d»T^) (1) 



a r 



Z A (1 + b A am q ) ^A^ + acA,d v 



P 



(2) 



2. DECAY CONSTANTS 



The results for the decay constants can be 
found in ||, as well as a discussion of the im- 
provement coefficients, but we present the main 
results again. The definition of the vector decay 
constants differs from the pseudoscalar and so we 
give it here: 



Ml 



(3) 



The decay constants were extracted from appro- 
priate ratios of correlation functions, and then 
extrapolated linearly in the light quark mass to 
K n and k s . The extrapolation to the B meson is 
motivated by HQS, of the form 



.i> il .u,iEri„i/,i/,v.U; ! ~ i§: -i§2 ) 



where i is P or Vfl and the function C(Mj) con- 
tains the leading logarithmic a s corrections to the 
scaling relation 

The results for the decay constants are shown 
in table 0. The large value of the decay constants 
is in part due to the choice of quantity used to set 

1 In this case /j = My / fy 
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the scale, fg varies from 218 MeV using 7- to set 
the scale to 186 using m p . This dominates the 
systematic uncertainty. The remaining system- 
atic error comes from difference between (3 = 6.2 
and (3 = 6.0, and choice of functional form for the 
heavy extrapolation (quadratic vs linear). Ac- 
cording to HQS, the $ functions for the P and V 
decay constants should agree in the heavy quark 
limit. However, the (3 = 6.0 data fail to sat- 
isfy this constraint, which may be indicative of 
large discretisation errors on the coarser lattice. 
To examine this further we considered, the KLM 
norm Q^Jand the dispersion relation. At {3 = 6.2, 
the KLM factor was close to unity and the disper- 
sion relation showed good continuum behaviour. 
However, at j3 = 6.0, the KLM factor has a fi- 
nite effect, and the dispersion relation shows a 
deviation from the continuum version. These all 
suggest that the (3 = 6.0 data are being affected 
by 0((arriQ) 2 ) errors. 

Table 2 



The decay constants. The first error is statistical, 
the second is systematic. 





P 


V 


B 
D 

B s 
D s 


218(5)+^ MeV 
220(3)± 2 4 MeV 
242(4) l*f MeV 
241(2)+ 3 ^ MeV 


22.6(0.7)+." 

7.5(0.1)±S;| 
20.9(0.4)+^ 

7.3(0.1)lg5 


B/B s 
D/D s 


1.11(0.01)™ 
1.09(0.01)18$ 


0.92(0.01)™ 
0.98(0.01)1°^ 



3. RENORMALISATION 

The matrix element of a P — > P transition can 
be parameterised in terms of two form factors, 

zf{pm\v»\p(pi)) Xatt (5) 

= f+(q 2 )(Pi +P2- A m 2<?) M + f (q 2 )A M 2q^ 

where q = pi — -pi and A m 2 is the difference of 
the square of the masses normalised by q 2 . Recall- 
ing that the vector current is conserved, we know 
that /+(0) must be one and therefore consider 
the forward, degenerate matrix element. Here, 
can then be measured directly by taking the 

2 modified to account for 0(a) improvement 



appropriate ratio of two- to three-point correla- 
tion functions. This can then be compared to 
Zy (1 + byamq) with the values of Zy and by de- 
termined by |5j, at the values of quark mass for 
which we measure Zy^ . This is shown in figure |l|. 
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Figure 1. The mass dependence of Zy. The 
vertical line shows the charm quark mass as de- 
termined from the heavy-light pseudoscalar. 

These data have already been shown for (3 = 
6.2 in §]. The ALPHA determinations of Zy 
and by used chiral ward identities, and so used 
zero or very small quark mass. Our results show 
excellent agreement with the mass-independent 
scheme even at large quark mass. 

4. THE SOFT PION RELATION 

HQS can be used to relate the scalar form fac- 
tor, /o to the ratio of the B to pion decay con- 
stant, 

/o(<?max) = -f- (6) 
J n 

A recent review J7J found that the majority of 
lattice calculations found this relation violated. 
However, the key issue for the SPR on the lattice 
is the chiral extrapolation. This is complicated 
by the dependence of the form factors on q 2 , it- 
self heavily dependent on the light quark mass. In 
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this calculation the active light quark (A), which 
occurs in the heavy-light current and the spec- 
tator light quark (S) which takes no part in the 
weak decay have different mass. Hence, 

A Taylor expansion of q 2 shows q 2 ~ M 

we know from PCAC that (Mp« 
the average quark mass fh — {tua 
the light quark mass dependence of the form fac- 
tor is 



flight \ 2 



(7) 

pf\ and 
fh where 
m s )/2. Thus 



(8) 



f(q 2 , m, n%s) = a + /3m 1//2 + 7m + Sins 

The second term in this expression makes the 
extrapolation difficult to control, as it is non- 
analytic and dm 1 / 2 /dm blows up as we approach 
the chiral limit, in the region where we have no 
data. 

Instead of extrapolating at fixed pion three- 
momentum we separate the q and chiral be- 
haviour of the form factor. We first interpolate 
in q 2 at fixed quark mass so that at each quark 
mass we have the form factors at the same mo- 
mentum value, and then extrapolate at fixed q 2 
without the troublesome second term. However, 
because we extrapolate at fixed q 2 , we no longer 
have data at g^ lax for physical quark masses. This 
method has already been successfully employed to 
calculate the differential decay rate for B — > it 
at j3 = 6.2 in Q. In this analysis, an extra- 
momentum channel is included (1 — > \/2||), and 
the range of q 2 altered to allow interpolation at 
fixed quark mass for both values of (3, in the same 
range of q 2 . 

The form of interpolating function is motivated 
by pole-dominance models. 

h 



(9) 



(l-q 2 /M 2 ) 

with kinematic constraint /o(0) = /+(0) imposed. 
We have also tried other functions for the inter- 
polation, including dipole/pole for /+ and fo re- 
spectively and with or without the kinematic con- 
straint. This is shown in figure |[ Also shown in 
figure |^ are other less well motivated functions, 
such as linear or quadratic in q 2 . As we are inter- 
polating any reasonable, smooth function should 
produce similar results. 
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Figure 2. Interpolation functions for heaviest k 
combinations at f3 = 6.0. The heavy quark has a 
mass around charm. 

For the extrapolation to the B meson mass, we 
are again guided by HQS, 
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where Ui = — 1, +1 when i = 
kinematic variable, 



v ■ k 



M 2 + M 2 - q 
2M 



, and v ■ k is the 



(11) 



The fixed q 2 method allows us to choose the q 2 
values at each heavy quark mass such that v ■ k 
remains constant during the heavy quark extrap- 
olation. 

We can now consider the momentum depen- 
dence of the form factors at the B scale. Com- 
bining HQS scaling relations and pole dominance 
models suggest the momentum dependence is a 
dipole for / + and a pole for /o (DPP). A more 
physical model has been suggested by Becirevic 
and Kaidalov (BK) @ 



Mi 2 ) 

fo(q 2 ) 



Cb(1 - a) 



c B (l - a) 
(l-q 2 /(3M 2 ,)- 



aq 2 /M 2 *) 



(12) 
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Figure 3. The momentum dependence of the form 
factors for B — > n. 

These are shown in figure ||. The value obtained 
for /o(<Zmax) can tnen be compared to the value of 
/b//tt- To try and estimate the systematic errors 
we considered the following; ro versus m p to set 
the scale, (3 = 6.2 versus — 6.0, different models 
for the fixed q 2 interpolation, quadratic vs linear 
heavy extrapolation. We see the SPR satisfied on 
the lattice, with large systematic errors, but with 
a tendency for / (g, 2 nax ) to lie below Jb/Itt- 

We can also check to see if the extrapolations 
in q 2 and heavy quark mass commute. In figure 
H the heavy extrapolation of Jb/ fn and /o(Smax) 
are shown after fo(q 2 ) has been extrapolated with 
a pole model to g^ ax . Again we see the SPR is 
satisfied. 

A recent calculation || using a NRQCD for- 
mulation of heavy quarks confirms the difficul- 
ties with the chiral extrapolations of the form 
factors and also adopts a fixed kinetic variable 
approach, which reduces their observed violation 
of the SPR. Another calculation using the 
relativistic heavy Wilson formulation, and a fixed 
kinematic variable approach see the SPR satisfied 
if a quadratic heavy extrapolation is used. 

The SPR is rather difficult to satisfy on the 
lattice because we have to extrapolate the light 
quark mass into a region where q 2 changes 
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Figure 4. The soft pion relation 

rapidly. This can be circumvented by using 
the fixed q 2 approach. However, lattice calcula- 
tions can reliably determine the differential decay 
rate for B — > tt away from g^ ax where the phase 
space is non-zero to extract V u h- This work 
was supported by EPSRC grant GR/K41663 
and PPARC grants PPA/G/S/1999/00022 
and PPA/GS/1997/00655. CMM acknowleges 
PPARC grant PPA/P/S/ 1998/00255. 
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